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ABSTRACT 

Previous treatments of ambipolar diffusion in star-forming molecular clouds do not 
consider the effects of fluctuations in the fluid fields about their mean values. This paper 
generalizes the ambipolar diffusion problem in molecular cloud layers to include such 
fluctuations. Because magnetic diffusion is a nonlinear process, fluctuations can lead 
to an enhancement of the ambipolar diffusion rate. In addition, the stochastic nature 
of the process makes the ambipolar diffusion time take on a distribution of different 
values. In this paper, we focus on the case of long wavelength fluctuations and find 
that the rate of ambipolar diffusion increases by a significant factor A ~ 1 — 10. The 
corresponding decrease in the magnetic diffusion time helps make ambipolar diffusion 
more consistent with observations. 

Subject headings: stars: formation - magnetohydrodynamics 

1. Introduction 



In the usual paradigm of low mass star formation, molecular cloud cores are supported by 
magnetic fields. In order for star formation to take place, the cores must lose magnetic support, 
and this loss of support is generally thought to take place through the action of ambipolar diffusion 
(Mouschovias 1976; Shu 1983; Nakano 1984; Shu, Adams, & Lizano 1987; Lizano & Shu 1989; 
Ciolek & Basu 2000, 2001). This general picture has support from observations, which suggest 
that ion- neutral drift does indeed occur in magnetized star-forming cores (e.g.. Greaves &: Holland 
1999). 

An important issue facing this standard scenario is the time scale required for magnetic sup- 
port to be removed from the cloud cores. As the observational picture comes into sharper focus, the 
number of observed cores without stars (e.g., Jijina, Myers, & Adams 1999) seems to be smaller than 
that predicted by most previous estimates from ambipolar diffusion (e.g., Ciolek & Mouschovias 
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1994, 1995; Lizano &; Shu 1989) by a factor of 3 - 10. In other words, loss of magnetic support by 

diffusion appears to be too slow, with a time scale a factor of 3 - 10 times longer than suggested by 
the observed statistics of cloud cores. However, these previous calculations neglected a dimension- 
less factor that depends on the mass to flux ratio of the cores (Ciolek &; Basu 2001). If the cloud 
cores have mass to flux ratios that approach the critical value, then the ambipolar diffusion time 
scale is significantly shorter than previous estimates. In particular, if the mass to flux ratio becomes 
supercritical, then the ambipolar diffusion time scale approaches zero. The need for this correction 
is bolstered by a recent compilation of Zeeman measurements of magnetic field strengths (Crutcher 
1999), which suggests that many cores may have mass to flux ratios near the supercritical value. 
This observed sample includes only 27 cores with relatively large masses; additional measurements 
are necessary to clarify the observational picture. 

In this study, we consider the effects of fluctuations on the mechanism of ambipolar diffusion 
described above. The time scale issue remains important and this work shows that ambipolar 
diffusion can operate more quickly in the presence of such fluctuations. In addition, because of the 
chaotic nature of the fluctuations, the ambipolar diffusion time scale will take on a full distribution 
of values for effectively "the same" initial states. 

Fluctuations are expected to be present in essentially all star forming regions. Molecular 
clouds are observed to have substantial non-thermal contributions to the observed molecular line- 
widths (e.g., Larson 1981; Myers, Ladd, &: Fuller 1991; Myers & Gammie 1999). These non- 
thermal motions are generally interpreted as arising from MHD turbulence (e.g., Arons & Max 
1975; Gammie & Ostriker 1996; for further evidence that the observed linewidths are magnetic in 
origin, see Mouschovias &; Psaltis 1995). Indeed, the size of these non-thermal motions, as indicated 
by the observed line- widths, are consistent with the magnitude of the Alfven speed (e.g., Myers & 
Goodman 1988; Crutcher 1998, 1999; McKee & Zwcibel 1995; Fatuzzo & Adams 1993). As a result, 
the fluctuations are often comparable in magnitude to the mean values of the fields (T. Troland, 
private communication) . 

Background fiuctuations can lead to a net change in the diffusion rate because magnetic dif- 
fusion is a nonlinear process. As many authors have derived previously (e.g., see the textbook 
treatment of Shu 1992), and as we present below, the (dimensionless) diffusion equation takes the 
schematic form 



where 5 is the magnetic field strength, /j, is the Lagrangian mass coordinate, and we have ignored 
density variations. Now suppose that the magnetic field fluctuates about its mean value on a time 
scale that is short compared to the diflFusion time (the time required for the mean value to change). 

We thus let b — ^ 6(1-|-^), where ^ is the relative fluctuation amplitude. In the simplest case in which 
the fluctuations are spatially independent, the right hand side of equation [1] is thus multiplied by 
a cubic factor (1 -|- ^)^. Although a linear correction would average out over time, this nonlinear 
term must always average to a value greater than unity and the corresponding diffusion time scale 
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grows shorter by the same factor. As an example, suppose the field spends half of its time at a 
value of twice its mean strength and the other half of its time near zero strength. Half of the time, 
the effective diffusion constant is thus larger by a factor 8, whereas the other half of the time, the 
diffusion constant is effectively zero. In this naive example, the mean diffusion constant is thus 4 
times larger due to the fluctuations. The goal of this paper is to derive a more rigorous argument 
for this time scale enhancement. 

This effect - changing diffusion time scales because of fluctuations - is well known in mathe- 
matical subfields. Simpler problems in which random noise flelds drive physical systems at different 
rates appear in a host of textbooks (e.g., Srinivasan k, Vasudevan 1971; Soong 1973). More re- 
cently, in the context of "stochastic ratchets", it has been shown that random fluctuations can 
drive a physical system to propagate "uphill", i.e., the opposite direction of its natural propagation 
in the absence of fluctuations (Doering, Horsthemke, & Riordan 1994). In astrophysics, stochastic 
aspects of magnetic field fluctuations have been considered in the context of cosmic ray propagation 
(e.g., Jokipii k. Parker 1969) and also in stellar atmospheres (e.g., Shore & Adelman 1976). Several 
previous papers have studied turbulent fluctuations in magnetically supported clouds, often by 
considering how the turbulence itself leads to field evolution in the absence of ambipolar diffusion 
(e.g., Kim 1997). The formation of cores through the dissipation of turbulence has been suggested 
(Myers & Lazarian 1998). Turbulence can also enhance the rate of ambipolar drift and may help 
explain the observed relationship between density and magnetic field strength, B oa p'^ (see Zweibel 
2001). In this current work, we consider ambipolar diffusion to be the main process that forms 
molecular cloud cores and study how fluctuations alter its rate. 

This paper is organized as follows. In §2, we reformulate the ambipolar diffusion calculation 
in a plane geometry, where we explicitly include fluctuations in both the magnetic field and the 
density field. We perform an analysis of the resulting set of equations in §3. We specialize to the 
limit of long wavelength fluctuations and apply the resulting formalism to astrophysical systems. 
We conclude, in §4, with a summary and discussion of our results. The case of short wavelength 
fluctuations is presented briefly in an Appendix. The most important outcome of this study is to 
demonstrate that fluctuations can lead to more rapid diffusion of magnetic fields in star forming 
regions and that the diffusion time scale takes on a distribution of values (rather than a single time 
scale) . 



2. Ambipolar Diffusion in the Presence of Fluctuations 

In this section, we modify the standard ambipolar diffusion derivation to include the effects 
of fluctuations. We consider the simplest case of a planar layer of molecular cloud material. The 
magnetic field lines are parallel to the plane and all quantities depend only on the height z above 
the midplane, 

B = B{z)x. (2) 
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Wc note that a host of two-dimensional axisymmetric calculations have already been done (Lizano 

& Shu 1989; Mouschovias & Morton 1991; Ciolek & Mouschovias 1994; Basu & Mouschovias 1994; 
Basu 1997, 1998; Ciolek &: Basu 2000). Although more idealized, this one-dimensional calculation 
is useful because it allows for analytical results and isolates the fluctuation effects that we are trying 
to elucidate (see also the discussion at the end of this section) . 

As derived previously (e.g., Shu 1992), with this basic configuration the equations of motion 
take the form 

| + |(H=0, (3) 

du du 1 d B'^\ 

|| = -4vrGp, (5) 

and finally 

dB d . d f dB\ 

5^ + 5i^^") = 5il4;^a7J- 

To close the system of equations, the ion mass density pi must be specified. The ion population 
depends on the complex balance between the ionization rate of the neutrals (primarily through 
cosmic rays), the subsequent production of molecular and metal ions, and the ion-electron recom- 
bination rate in the presence of grains. A reasonable (standard) approximation for the ionic mass 
density in molecular cloud environments is given by 

Pi = Cp'/', (7) 

where C = 3x 10~^^ cm~^/^ g^^^- However, a new observational study (Caselli et al. 2001) indicates 
that the ionization fraction can be significantly smaller than suggested by equation [7] with this 
value of the constant C (see also Ciolek & Mouschovias 1998). This reduction in ionization also 
acts to speed up the ambipolar diffusion rate and helps alleviate the time scale problem of interest. 
Notice, however, that this effect is independent of our present work and can be easily incorporated 
by using a different value of the constant C. 

Next, we introduce the fluctuations through the ansatz 

B^B{l + and p^p{l + r)). (8) 

The relative fluctuations and r] obey distributions, which ultimately determine the effects of 
these variations. Although we will specialize to particular distributions later on, we now keep the 
analysis as general as possible. In any case, the fluctuations are assumed to have zero mean so that 
the quantities B and p appearing in the equations of motion can be thought of as time-averaged 
quantities and play the same role in this generalized calculation as they do in previous treatments. 

Following many previous treatments, we rewrite the problem in terms of Lagrangian coordi- 
nates so that the basic variable is the surface density of neutrals between the midplane {z = 0) and 
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a height z. This change of variables takes the form 

a = I p{z',t)dz'. 
Jo 



(9) 



Notice that this Lagrangian coordinate measures the distance in terms of the mean density p rather 
than the fuh fluctuating density field p{l + rj). In terms of this new coordinate, the equation of 
continuity becomes 

f)r 1 

(10) 



dz 



1 
p 



For an isothermal equation of state (P = a^p), the force equation can be written in the form 



1 d 
1 + 1] da 



Svr 
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Since ambipolar diffusion is generally much slower than gravitational collapse, the cloud is expected 
to evolve near a magnetohydro dynamic equilibrium state. The inertial term in the force equation 
can therefore be justifiably ignored, leaving an expression which can be integrated to yield the 
quasi-magnetohydrodynamic equilibrium condition. In this case, however, the integration can only 
be carried out after a suitable averaging of the fluctuations has been done (see below) . Finally, the 
magnetic field evolves according to a nonlinear diffusion equation of the form 



d_ 

at 



d 



p(l + r/)J l + ?7 5(7|47r7CpV2(l + r/)3/2a(7 



[5(1 + 0] 



(12) 



where the partial derivatives are taken with respect to constant t or constant a. 



Following the notation of Shu (1983, 1992), we introduce a dimensionless surface density p, 
volume density p, magnetic field 6, vertical coordinate y, and time r as follows: 



a = GooP , 
27rGc72 

P= 2^P^ 

a"' 

z = — y, 



27rGc7o 



and 



t = 



2(27rG)V2_ 



27rGc7o 



r . 



(13) 

(14) 
(15) 
(16) 

(17) 



The closed set of equations which describe ambipolar diffusion can now be written as a continuity 
equation 

dy_ 

dp p ' 



(18) 
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a force balance equation 



y Id 

p{l + r])dy' + 7— ^b(l + n)] + 
l + r]dfi 



1 d 



and a magnetic diffusion equation 



d 



dr [p{l + 77) 



1 + 7/ - 



6^(1 + 



0, 



(19) 



(20) 



These three equations, in conjuntion with the specification of the fluctuations, constitute the prob- 
lem to be solved. 

As is well known, this one-dimensional slab model is unrealistic in the sense that gravity 
"saturates". In other words, in Lagrangian coordinates, our model has a constant gravitational 
field strength for a given column density. The gravitational field depends on the mass shell but 
not on the distance from the midplane. Hence, a Lagrangian observer would find no increase in 
the gravitational field strength as the slab is compressed, unlike the more realistic case of higher 
dimensions. This point is important because ambipolar diffusion is driven by self-gravity and the 
slab model lies in the regime of "weak gravity" (Mouschovias 1982). However, the enhancement of 
the time scale that we study here arises from the nonlinearity of the equations and will be present 
in any geometry. To iUustrate this claim, wc derive an analogous formulation for a cylindrical 
geometry in Appendix A. In particular, we find that the enhancement factor (due to fluctuations) 
takes the same form as that derived above. Although the absolute value of the ambipolar diffusion 
time scale depends on the geometry, the enhancement factor calculated here does not. 



3. Analysis 

In this section, we find solutions to a limited version of the problem derived in the previous 
section. Unfortunately, a full solution to equations [18 - 20] is beyond the scope of this paper and 
several important simplifications are necessary to make further progress. 

In a complete theory, one would also derive equations of motion for the fluctuations, represented 
here by ^ and 77, and one would solve for their temporal and spatial dependences in a self-consistent 
manner. Such a calculation would require a theory of how MHD waves arc produced, and how 
they cascade into turbulence. Because we have no working a priori theory of turbulence, however, 
our first simplification is to assume viable fluctuation distributions rather than calculate them. In 
other words, we adopt a semi-empirical approach. These fluctuations are observed, and we can use 
observations to constrain their form, but we do not explicitly calculate their behavior. 

Regardless of what distributions are chosen for the fluctuations, another complication arises. 
The equations of motion derived above describe the evolution of one particular realization of the 
problem. The usual way to solve such a problem is to set up the initial state and step forward in 
time (either analytically or numerically) by sampling the values of the fluctuations ^ and rj from 
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their (presumably known) distributions. As time progresses, the field evolves and one obtains a 
solution. However, this approach - which we will adopt - does not provide a full solution to the 
problem. Because the fluctuations obey a distribution, a probabilistic description of the overall 
evolution problem is necessary. In particular, the time evolution of the magnetic field will have a 
distribution of possible solutions. A full solution to the entire problem (see Doering 1990) would 
calculate the full distribution of possible time behaviors for the magnetic field (and other fluid 
variables). Because this issue is compounded for the case of short wavelength fluctuations, we 
focus our analysis on long wavelength fluctuations. Wc briefly discuss the case of short wavelength 
fluctuations in Appendix B, which shows that additional simplifying assumptions are necessary to 
make progress. 



3.1. Time Averaging 

In order to describe the average behavior of the fluid fields, we introduce a method of inter- 
mediate time averaging. Previous calculations for slab models suggest that ambipolar diffusion 
occurs on a time scale tad ~ 10 — 15 (Nakano 1984, Shu 1983, Mouschovias 1983), although this 
time scale is significantly shorter when the mass to flux ratio of the cores approach supercritical 
values (Ciolek & Basu 2001). For comparison, observations suggest that the real value is closer 
to TAD ~ 1 — 3. The fluctuations themselves presumably occur within the MHD regime and thus 
occur on time scales Tmhd ~ ^/vA, where A is the length scale of the fluctuation and va is the 
Alfven speed. We can write the length scale in dimensionless form through the relation 

27rG(j 00"^' ^ 

where we expect x ~ 1 — 10 for the long wavelength limit considered below (see §3.3). The 
dimensionless time scale for which the fluctuations vary is thus given by 



Tmhd 



7C 



2(27rG)V2 



-1 



-X, (22) 

VA 



where the dimensionless parameter [7C/2(27rG)^/^] ^ 8 for our adopted values of C and 7. For 
this approach to be consistent, the fluctuation time scale must be much shorter than the ambipolar 
diffusion time. This requirement implies a constraint of the form x ^ 8{vA/a)TAD, where wc have 
used standard values for 7 and C (see also Zweibel 1988). Notice that the ratio va/ci is typically 5 
- 10, so this constraint is usually easy to satisfy (see, however, the discussion below). 

To find time-averaged quantities, we must average over an intermediate time scale tq that 
obeys the ordering 

Tmhd < To < TAD ■ (23) 

Again, we require Tmhd tad so that an intermediate range of time scales exists. Further, we let 
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brackets (...) denote time averaged quantities so that 

1 r'^ 

{p{l + rj)) = - p{l + 7])dT' ^p, (24) 

To Jo 



and 



(6(1 +0) = - r + ^)^^' ~ ^ • (25) 

To Jo 

In equating the time averaged quantities with p and b, we are ignoring errors of O [tq/tad]- For 
the reasonable choice of taking the averaging time scale to be the geometric mean tq = ^/Tn^hdTAD, 
the relative error becomes O [(rmhd/TAD)^''^]- 



3.2. The Quasi-Equilibrium State 

Next, we need to consider the equilibrium states. The time scale over which fluctuations change 
(roughly given by the MHD crossing time) is much shorter than the collapse time scale (roughly 
given by the sound crossing time) for cloud cores undergoing ambipolar diffusion. Because the 
fluctuation time scale is also much shorter than the ambipolar diffusion time scale, the slab is 
expected to evolve in a quasi-static equilibrium state supported primarily by magnetic pressure. 
Since the force resulting from magnetic pressure is nonlinear, however, the fluctuations can play 
an important role in supporting this state. The quasi-equilibrium condition is found by integrating 
the time-averaged force equation [19] to obtain 

Kb^+p=l-H^ where K = ( ^ ^\ \ . (26) 

\ (1 + ??) / 

In this geometry (used here to make the problem more tractable), an equilibrium state exists even 
in the absence of magnetic fields and full gravitational collapse is not possible. While this geometry 
does not reflect realistic conditions in molecular cloud cores, it nicely illustrates the effects of 
fluctuations on the ambipolar diffusion problem. 

In obtaining the above form, we have assumed that the fluctuations are relatively well-behaved 
so that they obey constraints of the form 

^ ^"^^-0 and /(l + o|^\~0. (27) 



1 + r]dfi/ \ dfx 

In order for these constraints to hold, the derivatives of the fluctuations must average to zero (as 
expected) and the derivatives must not be correlated with the fluctuations themselves. In other 
words, the fluctuations must be both spatially and temporally symmetric. Ambipolar diffusion, as 
expressed by equation [20] , then evolves the quasi-equilibrium state from a magnetically supported 
configuration to a thermally supported one in the slab model. 
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3.3. Long Wavelength Fluctuations 

We now consider the effects of fluctuations on the ambipolar diffusion process. As shown below, 
these effects vary considerably with the spectrum of the fluctuations. To start, we assume that 
the fluctuations in density and the magnetic field are independent of each other. We also consider 
the simplest case in which the fluctuations have long wavelengths. In the limiting case where the 
fluctuation scale A is much larger than the typical length scale R in the cloud, the fluctuations can 
be considered to be spatially independent in the diffusion equation [20]. This regime of parameter 
space is defined by the constraint 

For the rest of this paper, we will specialize our analysis to long wavelength fluctuations that obey 
this constraint. 

Although nature can also support short wavelength fluctuations, a rigorous treatment of their 
effects is made difficult by several issues: (1) In stochastic differential equations, the solutions 
depend on the manner in which various limits are taken (sec Doering 1990). For short wavelength 
fluctuations, many different equivalent ways of taking the appropriate limits are possible and no 
unique solution exists without further specification. In this context, we would need to understand 
the origin of turbulent fluctuations to provide further specification. (2) In numerical treatments 
of the diffusion problem, the diffusion equation can be unstable on short length scales. Many 
numerical treatments allow for (small) errors on short size scales, but provide the correct global 
behavior (see Press et al. 1986). To adequately follow short wavelength fluctuations, one needs a 
numerical method that adequately resolves all spatial scales. In light of these difficulties, we focus 
on long wavelength perturbations in this paper. For purposes of illustration, however, we consider 
a representative approach for short wavelength fluctuations in Appendix B. 

In order for the fluctuations to affect the diffusion process, the effective diffusion coefficient 
must sample a range of possible fluctuations during the ambipolar diffusion time tad (see equation 
[17]). If the fluctuations are MHD in origin, their time scale is given roughly by the Alfven crossing 
time. Sufficiently rapid fluctuations thus impose the constraint 



"^rnhd 



2(27rG)V2 



—TADX'^ > 1 , (29) 
a 



where tad is the (non-dimensional) time interval required for diffusion to take place and where 
"^mhd is given by equation [22]. The quantity Np thus represents the number of times that the 
distributions are independently sampled during the course of the diffusion process. The time 
required for the field strength to decrease by a factor of ~ 2.7 is usually tad = 5—10. The ratio 
va/o. is typically 3 - 10, the length scale x = 1 — 10, and the dimensionless ratio in square brackets 
is about 8. For reasonable values of the parameters, the left hand side of the above inequality lies 
in the range Np = 20 - 800, a wide range but always comfortably greater than unity. Furthermore, 
we expect the variations from case to case to differ from the expectation values by relative sizes 



-10- 



— 1/2 

of order ~ Np ' due to incomplete sampling of the distribution functions during the diffusion 
process. This relative variation is thus expected to be ~ 0.04 — 0.2. 

When constraints [28] and [29] are satisfied, the fluctuations are independent of the spatial 
derivatives appearing in the diffusion equation [20]. Before time averaging the diffusion equation 
[20], we rewrite the equation in a specific form to simplify the averaging procedure. In particular, 
we write out the time derivatives, collect the terms on the right hand side, and multiply by one 
factor of (1 + r/) to obtain the form 

^p[^~''|TTi] " {l + ^f/''dii{p^^ ' ^^^^ 
where the dots represent time derivatives. With the diffusion equation written in this form, the 
time averaging procedure removes all fluctuations except for an effective diffusion constant on the 
right hand side, and the diffusion equation becomes 

Notice that this effective diffusion constant D would be unity in the absence of fluctuations. 



3.4. Results for particular distributions 

After time averaging, the diffusion problem is a rescaled version of the one solved previously. 
The quasi-static equilibrium state is given by equation [26], which contains the factor K due to 
the pressure provided by fluctuations. In the diffusion equation, the effective diffusion constant is 
larger by the factor D. If we rescale the magnetic field strength according to 6 = VKh, then the 
quasi-static equilibrium equation [26] takes its standard form. The diffusion equation also takes 
its standard form if we rescale the time coordinate according to r = {D /K)t. In other words, the 
ambipolar diffusion process speeds up by a factor A. = D/K which is determined by the distribution 
of the fluctuations. To make further progress, we thus have to specify the distributions. 

The magnetic field fiuctTiations ^ and the density fluctuations r] follow normalized distributions, 
/(^) and g{rj), with zero mean, i.e., 

/ fim = 1 , / f^d^ = , j giv)dv = 1 , and J g{i^)rjdi^ = , (32) 

where f{^)d$, is the probability that a fluctuation in the magnetic field has an amplitude between 
^ and ^ + (f^ (the function g{r]) is defined similarly for density fluctuations). For a given choice 
of distributions /(^) and g{ri), the expectation value of the diffusion constant D, the equilibrium 
factor K, and the time scale correction factor A = D/K can be calculated directly. 

For purposes of illustration, we first consider the fiuctuations to have uniform (flat) distribu- 
tions with amplitude A < 1. In other words, 

m = ^ for -A<^<A, (33) 
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and ^ 

9i^) = 2A ^^"^ -A<ri<A. (34) 
We are implicitly assuming that the density fluctuations are independent of the magnetic fleld 
fluctuations. With this particular choice for the distribution functions, the enhancement factor can 
be written in the form 

^ 6(1 + ^^)(1-A)-V^-(1 + A)-V^ 

3 + ln(l + yl) - ln(l - yl) ' ^ ' 

The result is plotted in Figure 1, which shows that the enhancement factor becomes substantial, A ~ 
5, as the fluctuation amplitude A approaches unity. In this particular treatment, the enhancement 
factor has a logarithmic divergence as — > 1. 

For comparison, if we ignore density fluctuations and consider only fluctuations in the magnetic 
field strength, the enhancement factor simplifies to the form 

where the first equality holds for (^^) = and where the second equality holds for a uniform 
distribution of fluctuations. For the benchmark case A = 1, the magnetic fleld fluctuations are 
comparable to the total field strength and we find a modest increase in the ambipolar diffusion 
rate, i.e., A = 1.5. In the extreme limit (i^^) oo, the enhancement factor A ^ 3 (for (^^) = 0). 
Thus, as expected, fluctuations in the magnetic fleld strength alone lead to smaller changes than 
when the density fleld has independent variations. 

We can also consider the case of density fluctuations that depend directly on the magnetic 
fleld fluctuations. For example, flux freezing arguments (Shu 1992) imply that B (x p'^, where k = 
2/3 for the usual case of spherical geometry and k = 1 for a one-dimensional cloud layer. If the 
flelds, including fluctuations about their mean values, obey such a flux freezing relation, then we 
have a correlation of the form {1 + rj) oc (1 + ^)^/'' and the enhancement factor becomes 

^_((i+o!:^ 

For our reference case of uniform distributions with amplitude A= l,we can evaluate this expres- 
sion to obtain A = (6k — 2)2/[(8k — 3)2^/^"^]. For the flux freezing exponent k = 1 appropriate for 
a one-dimensional layer, for example, we flnd A 1.13. 

When the fluctuations in magnetic field strength and density are correlated, the enhancement 
factor is generally smaller than in the absence of correlations. The special case of k = 1/2 is 
particularly interesting: As shown by equation [37], the value k = 1/2 leads to A = 1, i.e., no 
net enhancement of the ambipolar diffusion rate. In other words, for this particular correlation 
between the density and magnetic field fluctuations, the problem reduces to its old (non-fluctuating) 
form. This value k = 1/2 is found in three-dimensional calculations of magnetic cloud models 
(Mouschovias 1976) and is consistent with Zeeman measurements of the magnetic field strength 
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in several molecular clouds (Crutcher 1999; Basu 2000). However, both the theoretical calculation 

and observational relation correspond to the mean field and not the fluctuations; one interpretation 
is that the relation B oc p^/^ is more of an upper envelope than a scaling law (Zweibel 2001). In 
any case, the correlation of the magnetic field fluctuations with the density fluctuations thus needs 
to be further specified. 

In general, MHD disturbances exhibit a rich variety of different possible behaviors, including 
the various kinds of correlations between density and magnetic field fluctuations. In the linear 
regime, for example, pure Alfven waves exhibit magnetic field fiuctuations but have no density 
variations. Purely acoustic waves have density variations but no magnetic field fiuctuations. When 
the magnetic field is perpendicular to the direction of wave propagation, then fiuids can develop 
magnetosonic waves with correlated magnetic and density fluctuations; this configuration also sup- 
ports slow (nonpropagating) modes that compress the matter (increase density) and displace the 
magnetic field. In real molecular clouds, MHD disturbances are a complicated and non-linear su- 
perposition of many different types of waves and other motions (one should keep in mind that these 
different types of waves propagate at different speeds). As a result, many different possible distri- 
butions of fluctuations are allowed by the observations. Without further constraints, the parameter 
space of possible fluctuations allows arbitrarily large enhancements in the ambipolar diffusion rate. 
This claim is substantiated by equation [35], which shows that A increases without bound in the 
limit ^ — 1 for one particular choice of distribution. To obtain bounds on A, we thus need to 
impose additional constraints, as discussed below. 



3.5. Application to Molecular Cloud Cores 

We can scale the parameters of this theory to observed molecular clouds in order to specify 
the expected regime of our derived results. As mentioned in the Introduction, molecular clouds 
have significant non-thermal contributions to their line-widths (Larson 1981; Myers, Ladd, & Fuller 
1991) and the implied non-thermal motions are often interpreted as MHD turbulence. Indeed, the 
amplitude of these motions is consistent with the expected Alfven speed in these regions. 

To constrain the allowed range of our fluctuations, we consider the observed line-widths to 
arise from the fluctuating part of the Alfven speed in a time-averaged sense. In other words, we 
assume that 

(A*.H4>„.^(|;^\if)-(|;) (38) 

where we have subtracted off the non-fluctuating part so that the non-thermal contribution to 
the line-width vanishes in the absence of fluctuations. Evaluating this expression in terms of our 
formalism for the fluctuations, we find 

{Av)%T = 2a^ — [K-l\= 2aoa2 - 1] , (39) 
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where a is the sound speed, K measures the amphtude of the fluctuations as defined by equation 
[26], and oq is the initial ratio of magnetic to thermal pressure. To evaluate the ratio b'^/p, we have 
used our equilibrium state with no fluctuations as a reference state. Since the observed line-widths 
are comparable to the expected Alfven speeds, we expect (Av « 2aoa? and thus we expect the 
quantity {K — 1) to be close to unity. In other words, scaling our formulation to observed molecular 
clouds implies that K 

We need to apply this result to constrain the possible values of A. For fluctuations obeying 
uniform distributions, the parameter if = 2 for an amplitude A 0.886; the corresponding value of 
the enhancement factor is A 2.26. For the special case of no density fluctuations, the enhancement 
factor is given by equation [36]; the constraint K = 2 implies that (^^) = 1 and hence A = 2 (in 
the symmetric limit where (^'^) = 0). In general, the inclusion of density fluctuations makes the 
enhancement factor larger. For typical star forming regions, our net result is that fluctuations 
increase the ambipolar diffusion rate by a factor of 2. 



3.6. Comparison to Numerical Results 

We now compare our analytic results with numerical simulations. Instead of working exclu- 
sively with the expectation values {D and K) , we can use the same distributions for the fluctuations 
^ and rj and then sample the distributions as we numerically integrate the equations of motion. 
Specifically, we use a time-averaged quasi-static equilibrium state (equation [26]) and numerically 
integrate the difi'usion equation [31] by sampling the distribution functions for the fiuctuations to 
specify the effective diffusion constant. Our basic numerical scheme is similar to that described 
in Shu (1983), although the time resolution must be significantly higher (and highly variable) to 
properly include the wide range of values for the effective diffusion constant. Each numerical simu- 
lation results in one particular realization of the time evolution. If wc perform enough realizations 
of the problem, the average time scale should be the same as that calculated from the expectation 
values above. 

To numerically follow the evolution of the cloud, we need to specify the starting condition. 
Following previous authors, we adopt the following modified standard family of initial states 
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where the ratio a^/K represents the initial magnetic to thermal pressure ratio, i.e., 
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Through the action of ambipolar diffusion, the fluid evolves toward a final state described by the 
above equations in the limit ao — ^ 0. 

In addition to recovering the expectation values for the diffusion time, the numerical treatment 
also finds the deviations about the mean value (the expectation value). These deviations depend on 
how often we sample the distributions to evaluate the effective diffusion constant. If we sample the 
distributions often enough, say at every numerical time step, then the numerically calculated time 
scale converges to the expectation values found earlier. As outlined in §3.3, however, we expect 
the fiuctuations to vary on the MHD crossing time scale (equation [22]). In dimensionless time 
units, we thus expect the diffusion constant to take on independent new values for time intervals 
longer than tx ~ O.la/vA ~ 0.02. For the sake of definiteness, we take tx = 0.02 for all of the 
simulations presented here. As we show below, however, the value of tx determines the width of 
the distribution of possible time scales; we can estimate this width analytically and then scale our 
results to other choices of tx- 

We compare our numerical results with the analytical formulation in Figure 2. Wc have used 
uniform distributions of the fluctuations with amplitude A (see equations [33 - 34]) and set ao = 10 
so that the magnetic pressure is ten times larger than the thermal pressure at the start of the 
calculation. In these simulations, the diffusion constant Dj = (1 +^)^/(l + r])^^^ is changed every 
time interval tx = 0.02 by sampling the distributions of and r]. The resulting time scale Tc for 
the magnetic field strength to decrease by one e-folding is plotted versus the fluctuation amplitude 
in Figure 2. The solid curve shows the expectation value for the time scale Te = Teo/A, where Teo 
is the time scale in the absence of fluctuations. The symbols show different realizations found by 
numerically integrating the problem. For this set of simulations, we have selected the amplitude A 
randomly as well. As expected, the numerical integrations agree with the analytic predictions of 
§3.3. In particular, the mean value of the distribution of time scales from the numerical experiments 
closely follows the expectation value. 

The width of the distribution (for a given amplitude A) increases with the amplitude of the 
fluctuations. As as result, fluctuations not only force ambipolar diffusion to take place more rapidly, 
on average, but they also allow the process to sample a range of time scales. This range of possible 
time behaviors, as illustrated in Figure 2, is determined (in part) by the time scale tx that specifies 
how often the effective diffusion constant takes on different values. In the limit tx 0, the 
ambipolar diffusion time scale approaches the expectation value and the width of the time scale 
distribution shrinks to zero. In the opposite limit, the effective diffusion constant takes on a single 
(but randomly selected) value for the whole evolutionary time. In dimensionless units, the c-folding 
time for ambipolar diffusion is = 5 - 10 in this regime of parameter space. With tx = 0.02, near 
the low end of the range estimated from MHD considerations, the effective diffusion constant takes 
on N = 250 - 500 different values during an evolutionary run. 

We can quantify the width of the distribution of time scales for a given distribution of fluc- 
tuations and a given amplitude. For purposes of illustration, we use uniform distributions of 



-15- 



fiuctuations with amplitude A = 0.886, which corresponds to the case in which the fluctuations are 
large enough to account for the observed non-thermal line widths in molecular clouds (see §3.5). 
Figure 3 shows the distribution of ambipolar diffusion time scales, calculated both numerically and 
analytically. The numerical results were obtained by running the aforementioned code 10,000 times 
to build up the distribution shown by the solid curve in Figure 3. For comparison, the dashed 
curve shows a gaussian distribution whose width is calculated as described below. The solid spike 
at r 12.5 represents the delta-function distribution that applies in the absence of fluctuations. 

The distribution of time scales r can be estimated analytically. The peak of the distribution 
(denoted here as (r)) is determined by the enhancement factor A calculated in §3.4. For the case 
shown in Figure 3, A = 2.25 and hence (r) = 5.54. We can also find the width of the distribution. 
For a given realization of the problem, the time scale r is determined by the effective value D of 
the diffusion constant, i.e.. 



where Np is the number of independent samplings of the diffusion constant (equation [29]), D is 
the expectation value, and the Dj are the particular choices taken during a given run. (Notice that 
in the limit Np — oo, we must have D ^ D.) For simplicity, we consider Np = {t)/tx to have 
a fixed value, although the exact value also varies from case to case (this complication is a higher 
order effect). The case-to-case variation AD = D — D \n the diffusion constant is thus proportional 
to the sum of a large number (here, Np ~ 300) of random variables C,j = Dj/D — 1. These variables 
Cj are constructed to have zero mean and follow a distribution given by the convolution of /(^) and 
g{rj)] we denote the variance of this Cj distribution as a. Because of the central limit theorem, the 
distribution of the composite variable AD takes a gaussian form. Similarly, distribution of time 
scales takes a gaussian form and its width is given by 



We can compare this formula with our numerical results for the particular case shown in Figure 
3 (uniform fiuctuations with A = 0.886). The dashed curve shows a gaussian distribution with 
this predicted width; the distribution is in good agreement with the histogram of results (the solid 
curve) from the numerical simulations. More precisely, we find that the expectation value for the 
ambipolar diffusion time scale Te = Teo/A ~ 5.54 is in good agreement with that found numerically, 
i.e., (r) ~ 5.62. The chosen distributions of fluctuations show that a ~ 2.51 and hence equation 
[45] implies that the distribution of time scales should have width {a) ~ 0.835. The numerical 
simulations imply almost the same value {{a) « 0.834). For all quantities that can be compared, 
the numerical results and the analytic predictions agree to within about one percent. 

As shown by equation [45], the width of the time scale distribution depends on the time scale 
Tx over which the fluctuations occur, or, equivalently, the number Np of independent fluctuation 
samples. The value chosen for our numerical simulations, tx = 0.02, is near the low end of the 
expected range. For comparison, the dotted curve in Figure 3 also shows the distribution for a 
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Tx = 0.055, which corresponds to Np = 100. Notice that as the fluctuation time scale tx grows 
even larger (and Np decreases), the distribution of time scales grows wider and eventually departs 
from a gaussian form. 

4. Discussion 

In this paper, we have explored how fluctuations in the background flclds affect the rate of 
ambipolar diffusion. These fluctuations force the magnetic fleld strength and the density to sample 
a distribution of values, rather than take on a single value at a given point in space and time. We 
have used a one-dimensional molecular cloud layer as a test problem to study the effects of such 
fluctuations. 

The flrst principal result of this paper is that the time scale for ambipolar diffusion is altered by 
these fluctuations. In particular, fluctuations drive ambipolar diffusion to take place more rapidly, 
with the time scale shorter by a factor A ~ 1 — 10. For typical conditions in molecular clouds 
cores, the enhancement factor is near the lower end of this range, A ~ 2 — 3, but much larger 
enhancements remain possible. The ambipolar diffusion time scale depends on the distribution of 
fluctuations. For the case of uniform distributions, for example, the ambipolar diffusion time scale 
varies with the amplitude A as shown in Figure 1. In general, the time scale also depends on the 
shape of the distributions. 

For a given distribution of fluctuations, the ambipolar diffusion time scale also varies from 
realization to realization (see Figure 2). Descriptions of the ambipolar diffusion process thus face 
an interesting complication, which is our second principal result: The time scale for loss of magnetic 
support takes on a distribution of values instead of a single value. Consider two identical molecular 
cloud regions and suppose they are laced with fluctuations following a given distribution. Because 
the two regions experience incomplete (and different) samplings of the fluctuation distributions 
as their magnetic fields diffuse outwards, they will not exhibit the same diffusion time. This 
feature is more general than its manifestation in this particular test problem. When a physical 
system contains an effectively random element - in this context through chaos and turbulence - the 
outcomes must be described in terms of a probability distribution. For our test problem (ambipolar 
diffusion in a cloud layer), the single value of the e- folding time Te is replaced by a distribution of 
values (see Figure 3). Furthermore, the distribution of possible time scales approaches a gaussian 
form; the most likely value for the time scale is shorter than the case without fluctuations by the 
enhancement factor K = DjK (see equations [26, 31]) and the width of the distribution is given by 
equation [45]. 

This effect on the ambipolar diffusion time scale has important implications for star formation 
in molecular clouds. These clouds appear to be supported by magnetic flelds and the observed 
magnetic field strengths are commensurate with this view. However, statistics of molecular cloud 
cores (with and without young stellar objects) argues that the (uncorrected) time scale for ambipolar 
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diffusion may be too long to account for the observations. This work shows that magnetic fields 

can diff^use more rapidly than previous estimates suggest. This speed-up, along with any other 
enhancements (e.g., Ciolek &: Basu 2001), can help account for the observed statistics of molecular 
cloud cores. Another complicating issue arises: Because the ambipolar diffusion time scale takes 
on a distribution of values, and this distribution can be rather wide if the fluctuations change on 
long time scales tx , some core regions will experience much faster diffusion rates than others even 
if they have "the same" starting conditions. In this regime of diffusion activity, the cores that 
actually form stars are those which evolve on the "fast" side of the distribution, whereas the cores 
that happen to live on the "slow" side of the distribution will fail to form new stars. 

This preliminary treatment of fluctuations, including their effects on ambipolar diffusion and 
star formation, remains incomplete in several respects. In this paper, we have separated the cal- 
culation of the diffusion process from the determination of the fluctuations. In particular, we have 
assumed a priori forms for the fluctuations to study their implications. In a complete treatment, 
one should calculate the fluctuations and their effects in a self-consistent manner. In addition, 
we have focused on long wavelength fluctuations and have not considered spatial gradients in the 
fluctuating part of the fields. Magnetic turbulence cascades down to small scales, however, so it is 
possible that fields fiuctuate at length scales smaller than our MHD condition. This complication 
should also be considered in future work. Our present treatment is limited to one-dimensional 
slab models so that magnetic tension is not included; two-dimensional simulations should be done 
in the future. Another classical problem is the heating of molecular cloud regions by ambipolar 
diffusion; the effects of fluctuations on this mechanism should be considered. Finally, the act of star 
formation provides a source of new turbulence, which drives new fluctuations and can affect the 
ambipolar diffusion rates of neighboring cores; this feedback effect should also be studied. In any 
case, fluctuations in both the magnetic and density fields introduce an effectively random element 
into the ambipolar diffusion process, and thereby provide a rich class of new behavior for further 
study. 
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APPENDIX A: Cylindrical Geometry 

In this Appendix, we consider the effects of fluctuations on ambipolar diffusion in a molecular 
cloud filament. We consider only the simplest case of magnetic field lines that are aligned with the 
axis of the filament and depend only on the radial coordinate r, i.e., we have 



B = B{r) z . 



(Al) 
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With this basic configuration, the equations of motion take the form 
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where we have defined u to be the radial (and only nonvanishing) component of the velocity, and 
we have made use of the relationship defined by equation [7]. In practice, these filaments will 
be subject to clumping instabilities in both the linear (Gehman, Adams, & Watkins 1996) and 
nonlinear regimes (Adams, Fatuzzo, k, Watkins 1994); in this derivation, however, we neglect this 
issue and focus on the effects of the cylindrical geometry. 

Next, we introduce the fiuctuations through the ansatz given by equation [8] and rewrite the 
problem in terms of a Lagrangian description of the dynamics (e.g., see §2). For the cylindrical 
geometry considered here, the relevant Lagrangian coordinate is the mass per unit length a along 
the filament within a radius r of the central axis, i.e.. 



a = I p{r',t)r' dr' . 
Jo 



(A6) 



Notice that the variable a differs from the true mass per unit length by a factor of 27r which has 
been omitted for simplicity. The original problem in the variables (r, t) is now transformed to one 
in new variables {a, t) and the derivatives transform according to 
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With this transformation, the equation of continuity becomes 
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For an isothermal equation of state, the force equation can be written in the form 
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and the nonlinear diffusion equation for the magnetic field becomes 



d 



^(1 + 

dt [p{i + v) 



1 d f r2s2(l + ^)2 d ^ 



Equations [AlO - A12] exhibit exactly the same form as their slab counterparts (equations [10 - 
12]). As such, the effects of fluctuations on ambipolar diffusion in a cylindrical filament will scale 
exactly as for the slab geometry considered in the main text. 

APPENDIX B: Short Wavelength Fluctuations 

We consider the effects of short wavelength fluctuations in this Appendix. As noted in the 

main body of the text, this analysis is complicated by the fact that the solutions to stochastic 
differential equations depend on the manner in which various limits are taken (Doering 1990). The 
formulation presented below thus represents one possible approach to the general problem. 

The most likely source of short wavelength fluctuations is MHD turbulence, which is present in 
most regions of molecular clouds. The MHD condition, already built into the ambipolar diffusion 
equations, requires that the neutral fluid remain coupled to the ions and to the magnetic field. 
Physically, this condition is met if the ion-neutral collision frequency fin = ^pi exceeds the frequency 
associated with the MHD turbulence. The latter frequency can be approximated by fmhd ~ "^a/A, 
where va = B /(iTrp)^/"^ is the Alfven wave speed and A is the length scale of the fluctuations. As 
a result, the coupling condition requires that 



X > 0.09- 
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where % is the dimensionless turbulence length scale as defined by equation [21]. 



If we assume that the fluctuations are both spatially and temporally symmetric (which means 
that ^ and rf are not correlated with their first order derivatives), then the following relations hold: 
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for all well-behaved functions F{^, rf). Under these conditions, the quasi-equilibrium state described 
in S3. 2 remains valid. 



For short wavelength fluctuations, expanding the diffusion equation using the same approach 
as presented in §3.3 yields the form 
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Wc note that the second and third terms on the right hand side vanish when they are time averaged 
because they take the form given by equation [4]. With this simphfication, the time-averaged 
diffusion equation reduces to the form 
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where D is defined in equation [31] and 
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For the case in which the fluctuations ^ and rj are uncorrected, the parameter G simpUfies to the 
form 

Similarly, for the case of perfectly correlated fluctuations, G simplifies to the form 
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The relative size of the two terms on the right hand side of equation [B5] ultimately determines 
the behavior of the magnetic field diffusion. Since the fluctuations ^ and 77 vary on a length scale 
X <C 1, whereas b and p vary on a Icngthscale 1 + ao ~ Kv^/a^ S> 1, a simple scaling analysis 
naively suggests that the second term (with coefficient G) would dominate over the first (with 
coefficient D). Upon closer inspection, however, we see that the derivatives of the fluctuations tend 
to cancel out, so that the relative sizes of D and G depend on the form of the fluctuations. In 
any case, however, this treatment does not yield an expression that can be scaled to the previous 
solutions with no fluctuations. 
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FIGURE CAPTIONS 

Figure 1. Enhancement factor as a function of fluctuation amplitude. The solid curve shows the 
factor A by which the ambipolar diffusion process is sped up by fluctuations, as a function of 
their amplitude A for uniform distributions. The dashed curves shows the enhancement factor 
for fluctuations in the magnetic field only, whereas the dotted curve shows the result for density 
fluctuations. Notice that fluctuations in the magnetic field produce a greater effect for relatively 
small amplitudes, but the density fluctuations are more important for larger amplitudes. 

Figure 2. Comparison of numerical and analytic results. The time scale Te is the e-folding time 
for the magnetic field strength to decrease. The time Te is plotted here as a function of fluctuation 
amplitude A for uniform distributions of magnetic and density fluctuations. The solid curve shows 
the expectation value of the c-folding time scale as calculated analytically in §3.3. The symbols show 
results from numerical integrations of the same problem for different samplings of the distributions 
(see §3.5), where we have used tx = 0.02 as the time scale over which the diffusion constant changes 
(see text). 

Figure 3. Distribution of ambipolar diffusion times for a cloud layer with uniform fluctuations 
of amplitude A = 0.886; this level of fluctuations is consistent with the non-thermal line widths 
observed in star forming regions. The solid curve (histogram) shows the result of 10,000 numerical 
simulations with different realizations of the fluctuations. The dashed curve shows the analytic 
prediction for the time scale distribution - a gaussian with a peak value given by the expectation 
value and with a width predicted by application of the central limit theorem. The dotted curve 
depicts a wider gaussian distribution that applies for longer fluctuation time scales (here, Np = 100 
or Tx ~ 0.055). In the absence of fluctuations, the cloud maintains a single value for its ambipolar 
diffusion time, as shown by the delta-function spike at Te ~ 12.5. 



1.2 



"I 1 1 1 1 1 1 1 1 1 1 1 — I — r 



1 - 



0.8 



.r i 

■ r 

// 
/ 



.■ 1} 



0.4 



0.2 - 



L 



/ 

L 







■J---" I 



/ 



tl 

\ \ 

\ 
\ 

\ 

in 



J L 








J I I I 



10 



